Quantum phase properties of two-mode Jaynes-Cummings model for 
Schrodinger-cat states: interference and entanglement 
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In this paper we investigate the quantum phase properties for the coherent superposition 
states (Schrodinger-cat states) for two-mode multiphoton Jaynes-Cummings model in the 
framework of the Pegg-Barnett formalism. We also demonstrate the behavior of the Wigner 
(W) function at the phase space origin. We obtain many interesting results such as there is a 
clear relationship between the revival-collapse phenomenon occurring in the atomic inversion 

a ■ 

^ ■ (as well as in the evolution of the W function) and the behavior of the phase distribution of 



both the single-mode and two-mode cases. Furthermore, we find that the phase variances of 
the single-mode case can exhibit revival-collapse phenomenon about the long-time behavior. 
We show that such behavior occurs for interaction time several times smaller than that of 
the single-mode Jaynes-Cummings model. 

PACS numbers: 42.50Dv,42.60.Gd 

I. INTRODUCTION 

One of the few exactly solvable fully quantum-mechanically models describing the interaction 
between the single-mode electromagnetic wave and two-level atom is the Jaynes-Cummings model 
JCM [1]. Various nonclassical effects have been reported for this system, in particular, when the 
field is initially prepared in the coherent light. For instance, the revival-collapse phenomenon 
(RCP) in the atomic inversion [jj, quadrature squeezing [3] and sub-Poissonian statistics [4]. Fur- 
thermore, JCM has been used as a source for generating nonclassical states in the framework of 
th e corona, _ -tt^fl.tt. more details related to JCM reader can consult the 
review paper, e.g., Q|. Actually, this simple model, i.e. single-mode interacts with two-level atom, 
has been extended to include multimode fields [?|, multilevel atoms Q], multiatom interactions [jj], 



and clamping [lij. On the other hand, several experiments have been performed lor (he JCM. 



e.g., isolated single atoms a single-mode two-photon Rydberg atom "micromaser" 
observation of RCP in the evolution of the atomic inversion [131 ■ 

n 



and 



The progress in the gen eration of quantum mechanical states, e.g. in a trapped ion nonde- 



molation measurement 



151 ]. conditional measurement 



, ..etc., encourages researchers to analyse 



the behavior of the JCM for these states. One of these states is the Schrodinger-cat state 
which can exhibit various of the nonclassical effects arising from the interference in phase space. 
These states can be represented (for the jth mode of convenience) as 



\ a )tj = N A\ a i) + e i\ ~ Q i)]> 

(1) 

= £ <#V>, 

n=0 

where \atj) is a coherent state of the jth mode with complex amplitudes ay; €j is a c-number takes 
on values, 0, 1, —1 denoting coherent states, even coherent states and odd coherent states. Also Nj 
is the normalization constant having the form 



N 2 = - (2) 

1 [l + e2 + 2e i exp(-2|a j | 2 )]' {) 

In the second line of (pQ) we have expressed the states as a linear combination of Fock states, where 
the probability distribution amplitudes are 

1 a n 

cW = iy i ex P (--|a i | 2 )^[l + (-lf6 i ]. (3) 



Throughout this paper we use the shorthand notation Cn] instead of Cn^(ej). The two-mode 
density matrix of the cat state ([T|), which includes 16 elements, can be briefly written as 

Pf(0) = \a) e \a) e ' £ '(a\ e (a\ 

(4) 

= PI + PM, 

where pM and pi are the statistical-mixture and interference parts whose exact forms are given in 
\y\ . As we mentioned above that for the single-mode cat states the origin of the occurrence of 
the nonclassical effects is in the quantum interference between the components of the states. For 
instance, coherent states are close to classical states, however, even (odd) coherent states exhibit 
significant squeezing (sub-Poissonian statistics) and their photon number distributions provide 
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oscillatory behavior. It is worth reminding that the evolution of the Schrodinger-cat states with 



the single-mode JCM has been treated in 



18 



20 



2l(. 



The phase properties of the JCM using the Pegg-Barnett formalism 



in, e.t 



23 



24, 



221 ] have been investigated 



251 ] , in particular, when the field is initially prepared in the coherent state. It 



has been shown that for the single-mode JCM the evolution of both the phase variance and the 
jhase distribution can carry certain information on the RCP of the corresponding atomic inversion 



23J. Moreover, the phase variances of the multiphoton single-mode JCM exhibit RCP in the 



course of the long-time interaction 



241 ]. Recently, the phase properties of the two-mode JCM 



(TMJCM) in the form of linear interaction rather than nonlinear version is studied in [2q]. In the 
present paper we investigate the phase properties of the TMJCM whose Hamiltonian is described 
by the multiphoton hyper-Raman process and the modes are initially prepared in the Schrodinger- 
cat states (|4|). We use Pegg-Barnett formalism to perform such investigation. There are various 
mechanisms controlling the evolution of the considered system, e.g., interference in phase space, 
entanglement between the two modes (boson-boson interaction) as well as between the two modes 
and the two-level atom. The aim of the present analysis is to investigate the influence of these 
mechansims on the behavior of the phase distribution. Needless to say that the dynamics of the 
system are more complicated than those of the single-mode JCM. Moreover, such investigation is 
motivated by several facts. For instance, entanglement is connecting with the quantum computers 
and quantum cryptography 2?|] ■ Additionally, decoherence is at the heart of quantum measurement 
theory 28|]. Also the phase distribution has been realised experimentally by means of the optical 
homodyne tomography j^J, which allows quantum phase mean values to be calculated via the 
measured field density matrix. 

In this paper we obtain many interesting results for TMJCM. For example, there is a good 
connection between the behavior of the atomic inversion and that of the corresponding phase 
distribution of both the single-mode and compound-mode cases. Also we show that under certain 
condition there is a relationship between the evolution of the atomic inversion and the W function 
at the phase space origin. Furthermore, the evolution of the phase variance of the single-mode 
case can exhibit RCP about the long-time behavior, which occurs on the time domain several 
times smaller than that of the single-mode JCM. We analyse the behavior of the TMJCM in the 
following order: In section 2 we give the basic equations and relations, which will be used in the 
paper. Section 3 is devoted to discuss the main results. Summary of the results is given in section 
4. 
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II. BASIC EQUATIONS AND RELATIONS 

In this section we describe the system and provide the main relations used in the paper. Such 
relations include the Hamiltonian of the system, the corresponding wave function and the calcula- 



22] 



tions associated with the Pegg-Barnett formalism for the system 

We start with the Hamiltonian, which describes the interaction of the two modes with the single 
two-level atom and can briefly be written as 



T = H f +H T , 



(5) 



where Hf and Hj are the free part and the interaction part, respectively. In the rotating wave 
approximation these quantities take the forms 



(6) 



Hj = <7(a|"' 1 (1^0+ + a^a^a- 
where a± and a z are the Pauli spin operators; ujj, (j = 1,2) and co a are the frequencies of cavity 
modes a,j and the atomic transitions, respectively; g is the atom-field coupling constant and kj 



are the transition parameters, i.e. they represent the number of photons involvec 
transition. The derivation of the Hamiltonian ([6]) from the first principle is given in 
([6]) generalizes various models given in the literature, e.g 



31 



in the atomic 



321 ] . and for k\ 



301 ] . Moreover, 
■ = 1 it can be 



interpreted as a cavity version of Raman scattering in which a mode 1 is the pump field, mode 



311 ] . In this case the system represents three- 



2 is the Stokes field and with no anti-Stokes field 
level atom in A configuration interacting with the two modes of the field under the assumption of 
the exact resonance. Applying adiabatic elimination of the second level the effective Hamiltonian 
reduces to the form of the usual JCM Hamiltonian but the single-mode field operator is replaced 
by the product of the two modes operators. 

We proceed by considering the exact resonance case ui a = k\<jj\ — k<i<jJi- The dynamical properties 
of the system can be obtained by means of the unitary operator of the interaction part as 



U I (T,0) = eM-^H I ) 



cos(TD) - i s ^p-H I: 

gD 



(7) 



where 



T = gt, 



D 



it* 



CI2 0,2 <7-|_(T_ -+- CL^ CL^ 0,2 "2 



(8) 
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It is worth mentioning that without the rotating wave approximation, i.e. involving the fastly 
rotating terms in the Hamiltonian, the fully quantum-mechanical system cannot be solved as in 
([7]) since the eigenstates of the Hamiltonian cannot be obtained in a closed form. In this case, 
numerical techniques such as path-integral approach [33| and perterbative approach [3j] have to 
be used. Inspite of these techniques are very sophisticated the obtained solution will be in an 
approximate form. Of course, in this case the results of the present paper are no longer valid. For 
instance, we can mention for the single-mode JCM without RWA the atomic inversion is phase 
dependent [34j]. 



Additionally, we assume that the two modes are initially prepared in the cat state (j4|) and the 
atom is in the superposition of the excited and ground states having the from 

\<p,(j>) = cos <p\+) + exp(zc/>) sin tp\— ), (9) 

where |+) and |— ) denote excited and ground atomic states, respectively, and <j), <p are the relative 
phases between |— ) and |+). Therefore the total initial state of the system is 



l*(O)> = |a) 6l ($|a> e3 <g)|p,0). (10) 
From (fTJl- pIJj) the dynamical state of the system is 



|¥(T)) = l7j(T,0)|tf(0)) 

oo 

= [Fi(n,m,T)\+,n,m + k 2 ) + F 2 (n,m,T)\-,n + ki,m)] , 

n,m=0 

where the time-dependent coefficients have the forms: 
Fi(n, m, T) = C n 

;m+k2 COS V COs(TA nm ) iC n+kum exp(i0) sin cp sin(TA 

n,m J > 

F 2 (n,rn,T) = C n+kum sinvJcos(TA nim ) exp(z^) - iC n ^ m + k2 cos ipsm(TA ntTn ), 

whereas 



(11) 



(12) 



A n , m = J im + k f: + h)l , (13) 
V n\m\ 

and C n ^ m = Cn^Cm , C^P are given by ([3]). From (TTTj) the entanglement between different com- 
ponents of the system is readily apparent. 
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The atomic inversion of the system can be expressed as 

= E { K m+k2 cos 2 <p - C* +kum sin 2 <p] cos(2TA n>w ) 

(14) 

+ sin sin(2v?)C nim+fc2 C n+fclim sin(2TA„ >m ) | . 

To calculate the joint phase distribution for the field we find the reduced density operator for 
the radiation field pf(T) by tracing over the atomic states as 

p f (T) = Tr a \p(T)] 

OO r 

E Fi(n,m,T)F 1 *(n , ,m / ,T)|n,m + A: 2 )(n / ,m / + fc 2 | (15) 



n,m,n',rfi' 



+-F2(ra, nt, T)F£(n', m' , T)\n + fei, m)(n' + fci,m'| 



In the first line of (|15p the notation Tr a means that we evaluate the trace over the atomic states. 

Now we give briefly the relations of the Pegg-Barnett formalism [2] , which will be used through- 
out the paper, however, more details about the technique with various applications can be found 
in 



351 ] . This formalism is based on introducing a finite (s + l)-dimensional space ^ spanned by 
the number states |0), |1), |s). The expectation values of the different hermitian phase operators 
can be evaluated in the finite dimensional space Vl/ and at the final stage the limit s — > 00 is taken. 
The hermitian phase operator of the single-mode case is defined as 

s 

^e w |e m )(e m |, (16) 

m=0 

where the states |O m ) are eigenstates of the phase operator (fT6|) and they form a complete orthonor- 

mal basis of s + 1 states in ^. Moreover, these states are restricted to lie within a phase window 

00 

between Qq and 27r+Go- For the system described by the density matrix p = C m C*,\m){m'\ 

m,m'=0 

the continuum (, i.e. s tends to infinity) phase distribution is defined as 

p(e) = iim™^i(e m |p|e m ) 

(17) 

00 

= ^7 E C mC\ exp[i(m - m')8], 

m,m! =0 

where m has been replaced by the continuous phase variable G. Once the phase distribution P(0) 
is obtained, all the quantum-mechanical phase moments can be calculated as a classical integral 
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over G through the relation 

e 1 ' p(e,T)de, i' = i,2,... (is) 

-7T 

The phase variance is defined as 

((a4>) 2 ) = ($ 2 ) - (4>) 2 . (19) 



Actually, the generalization of the relations (116p - (|19j) to the two-mode version is straightforward. 
In this respect the joint phase distribution associated with the reduced density matrix (|15p can be 
evaluated as 

1 x 

P(©1, 2 , T) = — =■ > C(n, m, n, m', T) exp[i(n' — n)0i + i(m' — m)02], (20) 

where 

C(n, m, n', m', T) = Fx(n, to, T)F*(n', to', T) + F2(n, m, T)F£ (n', m 1 , T) 

(21) 

= [Cn,m+k 2 Cn',m'+k2 cos2 V 9 + Cn+ii.m^n'+t^m' V 5 ] c0S [P(A n ,m ~ A n / )jn /)]. 

For the sake of simplicity in (|2ip we restrict the calculation to the case <f> = . Also one can easily 
prove that 

r r p(e 1} e 2 , T)de lC ie 2 = 1. (22) 

J —TV J — IT 

The single-mode phase distribution of the jth mode can be obtained from (|20p via the relation 

r p(e r ,e j ,T)dQ f = p(e :i ,T), j^f. (23) 

1/ — 7T 

It is worth mentioning that (I20j) can be expressed in terms of the single-mode phase distribution 

as 



P(0i, 2 , T) = ^{2^ [P(6i, T) + P(0 2 , T)] - 1 



00 00 



+2 Yl Yl C + (n, m, n' , to', T) cos[(n — n')@i] cos[(to — m')0 2 ] (24) 

n>n' m>m' 

+C-(n, to, n', m',T) sin[(n — n')©i] sin[(m — m')©2] |, 

where 

C±(n, to, n', to', T) = C(n', to, n, to', T) ± C(n, m, n', to', T). (25) 
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Expression (|24p is relevant for numerical tasks. 

To investigate the fluctuation in the phase distribution for the system under consideration we 
have to evaluate different moments for the phase operators, which give 



oc 



^ -j^n/ +m' — m — n 



<<M> 2 > = - E C( ni m ) n'y i T)ig^ 
The sum- and difference-phase operators are defined as 



(27) 



where "+" and "-" denote sum-phase and difference-phase operators, respectively. Therefore, the 
sum- and difference-phase variances can be expressed in the following formula 



((A<l ± ) 2 ) = ((Adx) 2 ) + <(A<I 2 ) 2 > ± h lt 2, 



(28) 



/il J 2 = 2[($ 1 $ 2 )-($ 1 )($ 2 )]. 

Moreover, as the mean-photon number and the phase are conjugate quantities. It is reasonable to 
investigate the fluctuation in the mean-photon number for the single-mode and compound-mode 
cases ((Ahj) 2 ) and ((Ah±) 2 ), where hj = ajaj and n± = n\ ± n 2 . 

On the other hand, we address the relationship between the behaviors of the phase distribution 
and the W function at the phase space origin for the system under consideration. The reasons 
for this will be clear shortly. The joint Wigner function can be easily evaluated using a technique 
similar to that given in 



21( as 



W( X i,X2,T) 



cxp(-|xi| 2 -|X2| 2 ) 
^ 



E (-i 

n,m,n' ,m' 



in 1 '+m\ / n-n\ / m-m' 2 " +m ~ 2 " ~" 



Xl X 2 



(29) 



+(-l)^F 2 (n, m, T)F* 2 (n', m', T)yJ L — ^ (2| X i | 2 )L™r m ' (2| X2 | 2 )_ , 

where Xj = + Wj-,3 = 1,2 and L^(.) are associated Laguerre polynomials of order n. The 
function of the single-mode case can be obtained from (|29p by tracing over the other mode. Now 
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FIG. 1: The atomic inversion (a z (T)) versus the scaled time T for ip = and the modes are prepared 
initially in different types of the initial states with |aj | = 5. (a) (ei, e%, k\, fo) = (0,0,1,1) (curve A), 
(1, 1, 1, 1) (curve B); (b) (ei, £2, fci, fc 2 ) = (1, 0, 1, 1) (curve A), (1, 1, 2, 2) (curve B). The curves B are given 
for (a z (T)) + 2. 

the W functions associated with the first mode (Wi), second mode (W2) and compound modes 
(W) at the phase space origin, i.e. Xj = 0) can be expressed as 



W 1 (0,T) = - V (-l) n [|F 1 (n,m,T)| 2 + (-l) fcl |F 2 (n,m,T)| 5 



7T 



n,m=0 



(30) 



W 2 (0,T) = - V (-ir[(-l) fc »|Fi(n,m,T)| 2 + |F 2 (n,m,r)| i 

7T ^ — ' L 



n,m=0 



(31) 



^(°' T ) = ^2 E (-^[(-l^l^iCn.m.rjp + C-l^lFa^m.T)! 2 



n,m=0 



(32) 



Finally, as it is obvious that the system includes various parameters, which make the analysis to 
be difficult. Thus we restrict ourselves to the parameters, which can give significant results. Also 
in the text the statement standard JCM means k± = 0, fc 2 = 1> the optical cavity field and atom 
are initially prepared in coherent state and in the atomic excited state, respectively. 
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III. RESULTS AND DISCUSSIONS 

In this section we investigate the quantities, which have been calculated in section 2. Such 
investigation includes the atomic inversion, W function, phase distribution and and the phase and 
mean-photon variances. 

A. Atomic inversion and Wigner function 

As it is well known for the standard JCM that there is a direct connection between the phase 



distribution and (a z (T)) 23]. Additionally, quite recently it has been proved that there is a 



relation between (& Z (T)) and the evolution of the W function at the phase space ori gin 211 ] . 



This relationship is supported by the recent developments in photon counting experiment [36] and 



trapped ion technique 



381 ] . where the measurements have been focused on the origin of the W 
function in the phase space. The natural question is: what would be the situation for the TMJCM? 
In this part we answer this question by investigating the behavior of the (a z (T)) and showing its 
connection with the W function for the system under consideration. 

Firstly one of the most interesting phenomenon, which is representative to the single-mode JCM, 
is the RCP in the evolution of the (<r 2 (T)). Much attention has been focused on this phenomenon 
because it provides evidence of the granularity of the radiation field. It is worth mentioning that the 
observation of RCP has been performed via one-atom mazer [r]]. Moreover, schemes for measuring 



RCP via homodyne detection [39(], photon counting experiment and homodyne tomography [4C 
have been reported, too. For more details related to such phenomenon reader can consult, e.g. 
On the other hand, RCP associated with the TMJCM is very complicated compared to that one 
associated with the standard JCM where for the TMJCM the positions of the revivals in the time 
domain are independent of the intensities of the initial modes provided that (nj(0)) = fij >> 1. 



Such behavior has been explained partially in 



311 ] . Also in : 



it has been shown that there are 



secondary revivals coming from the complicated interferences, which are resulting from the double 
summations in the atomic inversion (cf. (|14p ). We should stress that the origin of the RCP in the 
(a z (T)) of the TMJCM is the strong entanglement between the bosonic and fermonic systems. To 
be more specific, restricting the discussion to the case k\ = &2 = 1, when the atom is classically 
treated the Hamiltonian ([5]) describes the up conversion process j^ll . 13 ] , which is representative 
well by switching the energy between the signal and idler modes. Conversely, when the field is 
classically treated, i.e. tij — > |ay| exp(— i(f)j), ([!]) reduces to two-level atom interacting with the 
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0, 0, 

FIG. 2: The phase distribution P(Q\) of the first mode when |aj| = 5, cp = 0, <f> = 0, k\ = fe 2 = 1 and for 
(ei,e 2 ) = (0, 0)(a) and (1, T)(b). For (a) T = 4.42 (solid curve-collapse time), 6.2999 (short-dashed curve- 
revival time) and 9.32 (long-dashed curve-secondary revival time). For (b) T — 1.8 (solid curve-collapse 
time), 3.22 (short-dashed curve-revival time) and 4.4 (long-dashed curve-secondary revival time). In (b) 
the star-centered curve is given for (ei, e 2 ) = (0, 1) and T = 6.3 (revival time). 

classical fields. In this case (a z (T)) exhibits steady-state Rabi frequency. In Figs. 1(a) and (b) 
we have plotted {a z (T)) against the scaled time T when the atom is in the excited state and the 
modes are initially prepared in different types of states as indicated. In Fig. 1(a) for curve A one 
can observe the revivals and secondary revivals representative to TMJCM. For curve B, where 
the interference in phase space starts to play a role, the revivals and secondary revivals are more 
pronounced than those shown in curve A and the behavior becomes more systematic, i.e. each 
revival is followed by a secondary revival. This can be explained as follows. For curve B the 
density matrix @ has two parts, which are pm and pi (cf. ((U)). Each of which provides its own 
RCP pattern and consequently they interfere with each other showing such behavior. Furthermore, 
it seems that the revival time for the curve A is two times greater than that for curve B. This can 
be partially explained for fij >> 1 and k\ = ki = 1. In this case the photon distribution for both 
initial coherent state and even coherent state possess Poissonian envelope. This means that the 
terms contributing effectively to the summations in (|14p are those for which the pairs (ai, CK2) and 
(n, m) are comparable. Now we assume that the arguments of cos(.) for initial coherent light is Ai 
with interaction time T and for initial even coherent light is A2 with interaction time T' , i.e. 

Ai = 2Ty/n(m + 1), A 2 = 2T'^2n{2m + 1). (33) 
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1.0 - 




FIG. 3: The phase distribution P(©i) of the single-mode JCM, i.e. k\ — 0, &2 = 1, when |cti | = 5,ip — 
0, = 0, ei = 1 and for T = 7 (solid curve-collapse time) and 16.4 (dashed curve- revival time). 

The revivals in (a z (T)) for both cases occur when Ai ~ 2m'ir, A2 ^ 2m' it, where m! is a positive 
integer. In the strong-intensity regime Ai and A2 can be simplified as 

Ai = 2TVfim~, A 2 = AT'Vfm. (34) 

Therefore, the relation between the revival times of the two cases is T' R = where the subscript 
R indicates that this time is associated with the revival pattern. On the other hand, in Fig. 1(b) 
the curve A is given for kj = 1 when one of the modes is initially in coherent states and the 
other is in even coherent states; while the curve B is given for kj = 2 and the modes are initially 
in even coherent states. The secondary revivals in curve A are more pronounced than those in 
Fig. 1(a). Nevertheless, the curve B, which is given for higher-order photon transition, exhibits 
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chaotic behavior resulting from the entanglement between different modes. This can be realised 
)y comparing this behavior with that of the atomic inversion of the two-photon-single-mode JCM 



211 ] . which exhibits compact and systematic revivals. From the above discussion we can conclude 
that the occurrence of the RCP and the secondary revivals in (a z (T)) for the TMJCM basically 
depends on the values of the transition parameters kj and on the type of the initial distribution of 
the modes. Nevertheless, it is insensitive to the values of the initial intensities of the modes. 

Now we draw the attention to the behavior of the W function, in particular, the evolution of 
the W function at the phase space origin. We start the discussion by assuming that the two modes 
are initially prepared in even coherent states and k\ = k 2 =odd number. In this case (130p -( [32]) 
reduce to 

W l (0,T) = -(a z (T)), W 2 (0,T) = —(a z (T)), W (0,T) = ^. (35) 

Expressions (|35p provide an important fact: The evolution of the phase-space origin of the single- 
mode W function can carry information on the atomic inversion of the system. However, the joint 
W function is localized in phase space giving maximum nonclassicality, i.e. maximum negativity. 
On the other hand, when the number of photons involved in the atomic transition is even, i.e. 
k\ = k 2 =even number, the W functions of the single mode and compound modes are localized in 
phase space: 

W l (0,T) = W 2 (0,T) = -, W(0,T) = -^. (36) 

Moreover, when k\ is odd, say, and k 2 is even we obtain 

W(0,T) = -Wx(0,T) = ^(a z (T)), W 2 (0,T) = -. (37) 

Finally, when k%, k 2 are odd and the two modes are initially prepared in a general quantum states, 
e.g. coherent states-restricting the discussion to the excited atomic state- ([30]) and ([32]) take the 
forms 

j oo 

W i(°> T ) = ^ E C tm+k 2 cos(2TA n , m + ror), (38) 

n,m=0 

W(0,T) =-Wi (0,0) W 2 (0,0), (39) 

where Wj (0, 0) is the initial value of the W function of the jth mode at the phase space origin. 
Expression (I38j) -apart from the prefactor l/w—is typical to that of the (a z (T)) but with additional 
factor, i.e. (—1)™, however, the joint W function is localized in the phase space. Actually, we have 
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FIG. 4: The joint phase distribution P(0i, 82) for |ay| = 5, = 0, (p — 0, fei = = 1, (ei, £2) = (1, 1) and 
for T = 1.8 (a), 4.4 (b) and 3.22 (c). 



numerically checked the behavior of f|38[) and found that it can provide typical curves to those of 
the corresponding atomic inversion. Therefore, for specific values of the interaction parameters 
the evolution of the W functions at the phase space origin represents the atomic inversion of the 
system, in particular, when k\ or ki (or both of them) is an odd number. This leads to the fact 
that for TMJCM there is a connection between the behavior of both the atomic inversion, the W 
function and the phase space distribution of the single mode as well as the compound modes, as 
we shall show shortly. This situation is similar to that of the standard JCM 231 ] . 



B. Phase distribution 



For the standard JCM it has been shown that there is a relationship between the time evolution 
of the atomic inversion, the phase distribution and Q function. For instance, the RCP is reflected 
in the behavior of the phase distribution of the radiation field [23j where the single-peak structure 
of the initial phase distribution splits into two peaks, which are rotated in opposite directions. At 
this moment (a z (T)) exhibits collapse pattern. Furthermore, as the interaction proceeds the two 
peaks of the distribution are collided showing single-peak form, however, (a z (T)) shows revivals. 
Here we show that similar behavior can occur for the single-mode and compound-mode cases of the 
TMJCM. Additionally, we investigate the influence of the higher-order photon transition kj > 1 
on such behavior. It is worth reminding that the behavior of the (a z (T)) of the TMJCM is rather 
complicated compared to that of the standard JCM (see Figs. 1). Actually, we found that the 
behavior of the first and second modes phase distribution are typical so that the discussion would 
be restricted to that of the first mode only. Furthermore, as it is difficult to treat the behavior of 
the phase distributions analytically we use the numerical technique (see Figs. 2-5). Figs. 2(a) and 
(b) are given for the single-mode phase distribution when the two modes are initially prepared in 
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0.60 □ 




FIG. 5: The phase distribution P(0i) of the first mode for T = 5, \atj\ — 5, (p = 0, <f> = 0, fci = hi = 2 and 
for (£1,62) = (0,0) (solid curve) and (1,1) (dashed curve). 

coherent and even coherent states, respectively, for collapse time (solid curve), secondary revival 
time (long-dashed curve) and revival time (short-dashed curve). In Fig. 2(b) the star-centered 
curve is given for the case (£1,62) = (0,1) and the corresponding revival time T = 6.3. In Fig. 
2(a) one can see that P(@x) exhibits two-peak structure for both collapse time (with two maxima 
around 0i ~ ±27r/3) and secondary-revival time (with two maxima around ®i ~ ±7r/2), however, 
those of the latter are border and shorter than those of the former. In fact such behavior (, i.e., 
the occurrence of the two-peak structure in the phase distribution) is quite similar to that of the 
even (odd) coherent states [2]. Furthermore, this behavior is representative to the standard JCM 
indicating generation of the cat states. Nevertheless, it is difficult to consider such conclusion here 
owing to the strong entanglement between different components of the system. For the revival time 
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P(@i) exhibits two wings around ©i — > ±7r (see short-dashed curve), i.e. it is incomplete peak. 
In Fig. 2(b), one can see that four-peak structure occurring through the collapse and secondary 
revival times for the initial even coherent states, which are combined at the revival time providing 
two-peak structure. Interestingly, when one of the modes is initially prepared in coherent state and 
the other is in even coherent state we have seen that for collapse and secondary revival times the 
behavior is typical to that of the two modes, which are initially prepared in the cat states (see solid 
and long-dashed curves in Fig. 2(b)), however, for the revival time we obtain single-peak around 
01 ~ and two wings as @i — > ±7r (see the star-centered curve in Fig. 2(b)). This indicates 
that the influence of the interference in phase space on the behavior of the phase distribution is 
more pronounced in the course of the revival time. Nevertheless, the influence of the entanglement 
on the behavior of the distribution can be realised by comparing Fig. 3, which is given for the 
single-mode JCM when the mode is initially prepared in the cat state, with the Fig. 2(b), in 
particular, the solid and short-dashed curves in the two figures. Such comparison shows that there 
is agreement and disagreement in the phase distribution of the two systems. For instance, the two 
systems exhibit four (two) peaks in the course of the collapse (revival) time, however, those of the 
single-mode JCM are narrower and higher than those of the TMJCM. The reason for occurring 
the same number of peaks in the distribution of the two systems is the term cos((n — n')@x). Also 
from (|35|) and (j38[) Figs. 2 can provide information on the evolution of the W function at the phase 
space origin. 

Now we turn the attention to the joint phase distribution. We have noted that the evolution 
of the P(Gi,02,T) reflects well the behavior of the (a z (T)) (see Figs. 4(a)-(c) for given values 
of the parameters, which are the same as those of Fig. 2b). Based on the fact P(©i,02,T) is 
symmetric in 0i and ©2 we have plotted Figs. 4 only over the region < 0j < vr. Generally, 
there is a similarity between the single-mode and compound-mode phase distribution. To be more 
specific, for initial even coherent light, P(0i, G2, T) exhibits two peaks over the specified region in 
the course of the collapse and secondary revival times, which reduce to single peak (recombination 
of the peaks) through the revival time. Moreover, the peaks, which are related to the secondary 
revival time, are broader than those of the collapse time (compare Fig. 4(a) and Fig. 4(b)). Also 
the influence of the interference in phase space on the behavior of the joint phase distribution is 
remarkable as multipeak structure. 

We conclude this part by shedding the light on the two facts, (i) For kj > 1 where the atomic 
inversion exhibiting chaotic behavior the phase distributions show multipeak structure, which is 
insensitive to the values of the interaction time (, i.e. steady state phase distribution). Of course 
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FIG. 6: The single-mode phase variance ((A$i) 2 ) of the first mode versus the scaled time T for \atj 
5,(p = 0,(j) = Q,k 1 = k2 = l, and for (ei,e 2 ) = (0,0) (a) and (1,1) (b). 



the number of peaks for the initial cat states are approximately two times greater than those for 
the initial coherent light (see Fig. 5 for given values of the parameters), (ii) As it is well known 
that for the standard JCM "coherent trapping" (, i.e. (a z (T)) ~ 0) can occur 44]. Similarly 
this can occur for the TMJCM under certain condition (cf. (|14p ). Nevertheless, this cannot be 
remarked in the behavior of the phase distribution as a result of the fact that the phase distribution 
is normalized (cf. (|22p ). Therefore, the relationship between the behavior of (a z (T)),P(Qj,T) and 
P(@i, ®2,T) occurs only when the two-level atom is either in the excited or ground state. Thus 
we have limited the discussion to these cases. Finally, we have found also that the behavior of 
P(0j,T) and P(@i,@2,T) are independent of the type of the initial cat states and the initial 
atomic states. 



Phase fluctuations 



For the standard JCM it has been shown that the phase variance exhibits collapses and revivals 
about the long-time interaction, and there appear main revivals and secondary revivals with differ- 



ent periods 



241 ]. The origin of this phenomenon is in the phase correlations between the different 



eigenstates of the photon number that make contributions to the orientated distribution of phase, 
and the contribution due to the correlations between the neighboring Fock states is much larger 



than the others 



24|. Actually, we have seen that this phenomenon can occur for the TMJCM for 
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interaction time shorter than that of the single-mode JCM (see Figs. 6-7 for given values of the 
interaction parameters). From Fig. 6(a) we observe that the RCP as well as secondary revival 
occur around 7r 2 /3, around the random phase distribution. Comparison of this figure with the 
Fig. 1 in [24| shows that the entanglement between the two modes in the interaction with the two- 
level atom makes the RCP in the phase variances occurring through interaction time several times 
smaller than that of the standard JCM. This behavior can be roughly explained in the following 
sense. We restrict the analysis to the case k\ = k\ = 1 when the two modes are initially prepared 
in the coherent states with strong intensities (, i.e. |oy| >> 1, j = 1,2) and the atom is in the 
excited state. From (I26h the phase variance of the first mode ((A<I>i) 2 ) takes the form 

n>n' m=0 

(40) 

x cospVm + l(Vn + 1 - y/ri + 1)]. 
It is worth mentioning that the expression of the phase variance of the standard JCM can be 
obtained from (|40p by dropping the summation Y^m=o an< ^ ^he notations related to the second 
mode. In this case the argument of the cos(.), i.e. the argument of the dynamical term of the 
standard JCM, is 

TZ n , n , = T(Vn~+l - vVTT). (41) 

Therefore, the phase variance can exhibit revivals, i.e. cos(.) provides maxima, only when 

TZ n ^ n i = In, I is integer (42) 

where T is the revival time, say, of the phase variance of the standard JCM. Now we draw the 
attention to the ((A<3E>i) 2 ) of the TMJCM given by (|40p . In the strong-intensity regime, i.e. x/fh = 
|ck2 | >> 1, we can apply the harmonic approximation technique [21] to evaluate the summation 
related to the second mode where {Cm) 2 has Poissonian distribution. In this case (I40p reduces to 

n>n' 

(43) 

x exp[-2msin 2 (^^)] cos^Ty/mZ^ +fhsm(^g-)}. 
Expression (|43p can provide revivals when the exponential term is maximum, i.e. 
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where T' denotes the revival time of the ((A$i) 2 ) of the TMJCM. Pr om (1421 ) and ( 1441 ) we arrive 
at 

r 



T. (45) 



4vm 

Expression (gSJ) shows that T' < T. In other words, the RCP in ((A<J>i) 2 ) for the TMJCM occurs 
for interaction time Ay^i times shorter than that of the standard JCM. We proceed that the RCP 
in ((A<E>i) 2 ) can occur through interaction time smaller than that shown in Fig. 6(a) by preparing 
the two modes initially in the cat states. This is shown in Fig. 6(b) and can be proved "roughly" 
using procedures as those given above. Moreover, the explanation used for the behavior of the 
atomic inversion of this case can be adopted here. 

On the other hand, for the evolution of the phase variances of the standard JCM there are two 
important cases, which are three- and four-photon transition, i.e. k\ = 0, hi = 3 and 4. The 
phase variances for these cases exhibit significant RCP, in particular, the four-photon transition 
case whose phase variance exhibits periodic revivals with period tt. Moreover, these revivals almost 
restore their initial amplitudes and widths, and they are independent of the intensities of the initial 
modes Surprisingly, we have noted-apart from the amplitudes of the revivals-that there is a 
one-to-one correspondence between the behavior of the phase variance of the standard JCM and 
that of the TMJCM. In other words, we have seen that the behavior of the single-mode phase 
variances of the cases ki + ki = 3 (provided that kj ^ 0, j = 1,2) and k\ = ki = 2 associated 
with the TMJCM are similar to those of the standard JCM. This information can be realised from 
Figs. 7(a) and (b). From Fig. 7(a) we can see a systematic behavior, i.e. each revival is followed 
by secondary revival and the initial revivals are restored. Furthermore, the revivals and secondary 
revivals are intensities independent (provided that |ay| >> 1) and they occur periodically with 
period tt/2, i.e. half of that of the standard JCM. The origin of such behavior is in the two- 
photon transition and in the interference in phase space (compare Fig. 7(a) with Fig. 4 in 24j]). 
Furthermore, the comparison of Fig. 7(b) with Fig. 3 in [24] is instructive, where one can realise 
the role of the entanglement and the interference in phase space. 

Generally, the occurrence of the RCP in the ((A$i) 2 ) cannot occur in the ((A&±) 2 ) where 
hi 2 7^ (cf. (I26j) ), We shed the light on the behavior of the mean-photon number variances. It is 
worth reminding that the mean-photon numbers and the phase operators are conjugate quantities 
(noncommuting operators). In Figs. 8(a), (b) and (c) we plot ((Ani) 2 )), ((An + ) 2 )) and ((An_) 2 ), 
respectively, for the given values of the interaction parameters. It is obvious that ((Ani) 2 )) and 
((An_) 2 ) exhibit RCP (see Fig. 8(a) and (c)) which is insensitive to the intensities of the initial 
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FIG. 7: The single-mode phase variance ((A$i) 2 ) of the first mode versus the scaled time T for \ctj\ — 
5,ej = 1,(^ = 0,0 = 0, (jfei.fe) = (2,2) (a) and (1,2) (b). 




FIG. 8: The mean-photon number variances versus the scaled time T for \aj\ = 5, ej = 1, <p = 0, <j> = 
0, (ki,k2) = (1, 1) and for the single-mode case (a), sum-mode case (b) and difference-mode case (c). 

fields. Moreover, this behavior is almost similar to that of the corresponding atomic inversion. In 
Fig. 8(b) we can see that the chaotic behavior is dominant. To be more specific, the behavior of 
((An + ) 2 ) is close to the steady state, ((An + (T)) 2 )) ~ ((An + (0)) 2 )). Comparison between Fig. 
8(b) and Fig. 8(c) shows that the correlation term, i.e. (nif^) — (ni)(n2), plays a constructive 
(destructive) role in the behavior of the ((An_) 2 ) (((An + ) 2 )). 

Now we can conclude that ((A<I>i) 2 ) can provide RCP, which is insensitive to the type of initial 
atomic state. We have to stress that the RCP in the atomic inversion and in the phase variances 
are completely different and cannot be compared. Finally, it is worth mentioning that RCP has 
been realised also for the phase variances of the two-photon down-conversion with quantized pump, 
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e.g. [45]. 



IV. CONCLUSIONS 



In this paper we have discussed the phase properties of the lossless multiphoton two-mode 
Jaynes-Cummings model for Schrodinger-cat states in the framework of Pegg-Barnett formalism. 
The investigation has included also the dynamical behavior of the W function at the phase space 
origin. We have shown that the evolution of the (& Z (T)) of the TMJCM is reflected in the behav- 
ior of the phase distribution where the splittings and overlappings of the phase distribution are 
accompanied by the collapses and revivals, respectively, in (a z (T)). This has been remarked in the 
behavior of both the single-mode and compound-mode phase distributions. Furthermore, in the 
course of the secondary revival in (a z (T)) the phase distribution exhibits behavior almost similar 
to that of the collapse case. Moreover, the nonclassical multipeak structure, which is representative 
to cat states, is remarked in the behavior of the phase distribution. Nevertheless, this behavior 
is insensitive to the type of both the initial cat state (, i.e. if the modes are initially prepared in 
even or odd coherent states) and the initial atomic states (, i.e. if it is in the excited or ground 
state). Additionally, we have shown that for higher-order photon transition, i.e. kj > 1, the phase 
distribution is in the steady state showing multipeak structure. Also we have proved that under 
certain condition there is a clear relationship between the atomic inversion and W function. Such 
relation leads to the fact: There is a relationship between the W(0, T) and the corresponding phase 
distribution. Finally, we have shown that the single-mode phase variances can exhibit RCP about 
the long-time interaction. 
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